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Many  high  temperature  oxides  and  low  temperature  polymeric  materials  transport  protons,  oxygen  ions, 
and  electrons  or  holes.  These  materials  are  candidates  as  membranes  for  fuel  cells  and  electrolyzers. 
This  manuscript  examines  non-equilibrium  steady  state  transport  through  such  membranes  under  the 
assumption  of  local  equilibrium.  A  simple  equivalent  circuit  analysis  is  given  for  transport  through  mixed 
proton,  oxygen  ion,  and  electron  (hole)  conducting  membranes.  The  cell  potential  can  be  described  in 
terms  of  transport  parameters  of  charged  species  (H+,  O2-  and  e  or  h)  and  internal  EMFs  given  in  terms 
of  chemical  potentials  of  neutral  species  {/i h2,  pLo2 ).  The  resulting  equation  for  cell  potential  is  similar  to 
the  Goldman-Hodgkin-Katz  (GHK)  equation  used  in  cell  physiology.  Transport  through  a  fuel  cell  based 
on  such  materials  is  examined.  Effects  of  electrolyte/electrode  interfaces  are  explicitly  included  in  the 
analysis.  Fluxes  of  H2  and  02  through  the  membrane  are  evaluated  at  open  circuit  and  under  load.  They 
obey  Onsager  reciprocity  relations,  inclusive  of  interface  effects.  The  analysis  also  shows  that  the  chemical 
potentials  of  H2  and  02  and  electric  potential,  tp,  exhibit  abrupt  changes  across  interfaces. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  aqueous  systems  of  interest  in  cell  physiology  and  other 
areas,  transport  of  electrically  charged  species  is  commonplace. 
Examples  include  transport  of  ionic  species  across  biological  or 
polymeric  membranes.  In  such  cases,  either  side  of  the  membrane 
has  aqueous  media  containing  ionic  species  in  solution.  Concen¬ 
trations  of  the  various  ionic  species  are  generally  different  on  the 
two  sides  of  the  cell  membrane.  The  cell  membrane  also  typically 
contains  water  (hydrated).  Transport  of  ionic  species  across  cell 
membranes  occurs  under  electrochemical  potential  gradients.  Elec¬ 
trically  charged  species  are  ionic  (cations  and  anions)  and  no  net 
electrical  current  transports  through  the  membrane.  Several  ionic 
species  can  be  in  solution  and  their  transport  occurs  to  varying 
degrees,  depending  upon  their  relative  concentrations  and  mobil¬ 
ities.  In  a  non-equilibrium  steady  state  with  concentrations  of  the 
various  species  in  the  two  reservoirs  assumed  fixed,  steady  state 
fluxes  of  ionic  species  are  established  across  the  membrane.  This 
leads  to  the  establishment  of  an  electrical  voltage  difference  across 
the  cell  membrane,  the  cell  potential,  which  can  be  described  in 
terms  of  the  concentrations  of  ionic  species  in  the  two  reservoirs 
and  their  transport  parameters  across  the  membrane.  The  resulting 
equation  is  known  as  the  Goldman-Hodgkin-Katz  (GHK)  equa¬ 
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tion  or  sometimes  as  Goldman  equation  [1,2].  The  GHK  equation 
assumes  equilibration  of  compositions  of  the  various  ionic  species 
at  liquid/membrane  interfaces. 

Aqueous  and  other  liquid  electrolyte  solutions  always  contain 
mobile  cations  and  mobile  anions.  Thus,  local  equilibrium  can  be 
described  in  terms  of  local  chemical  or  electrochemical  potentials 
of  the  ions.  For  example,  for  an  ionic  salt  of  the  type  MmXn  dissolved 
in  water,  the  dissociation  reaction  at  equilibrium  is 

MmXn  &  mMn+  +  nXm~  ( 1 ) 

where  the  valence  of  the  cation  is  n  and  that  of  the  anion  is  -m.  The 
corresponding  local  equilibrium  may  be  given  by 

MmXn  =  WlHMn+  +  n[AXm-  =  m/XMn+  +  njlXm  (2) 

where  the  /x2s  are  the  respective  chemical  potentials  and  /x2s  are 
the  electrochemical  potentials,  the  latter  given  by 

Ai^f+ZjF#  (3) 

where  z2-  is  the  valence,  F  is  the  Faraday  constant,1  and  0  is  the  local 
electrostatic  potential.  Although  not  commonly  used,  it  is  under¬ 
stood  that  local  equilibria  may  additionally  be  given  with  respect 
to  the  neutral  species  by  reactions  of  the  type 

M  &  Mn+  +  ner  (4) 


1  All  chemical  potentials  and  electrochemical  potentials  are  defined  on  a  molar 
basis. 
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Nomenclature 


List  of  symbols 

I  electrolyte  thickness  (cm) 

r ce  =  ^  cathode/electrolyte  interface  electronic  area  specific 
resistance  (£2  cm2) 

r “  =  %  anode/electrolyte  interface  electronic  area  specific 

°e 

resistance  (£2  cm2) 

rf  =  electrolyte  electronic  area  specific  resistance 

(£2  cm2) 

rc2  =  cathode/electrolyte  interface  oxygen  ion  area 

o2~ 

specific  resistance  (£2  cm2) 

ra2_  =  anode/electrolyte  interface  oxygen  ion  area  spe- 

0  °o2— 

cific  resistance  (£2  cm2) 

r®*2_  =  -Jj—  electrolyte  oxygen  ion  area  specific  resistance 
(£2  cm2) 

r£+  =  Jf-  cathode/electrolyte  interface  proton  area  specific 
resistance  (£2  cm2) 

r“+  =  anode/electrolyte  interface  proton  area  specific 
resistance  (£2  cm2) 

rf+  =  electrolyte  proton  area  specific  resistance  (£2  cm2) 

°H+ 

Jq 2—  oxygen  ion  current  density  through  the  membrane 

(A  cm-2) 

/H+  proton  current  density  through  the  membrane 

(A  cm-2) 

le  electronic  current  density  through  the  membrane 

(A  cm-2) 

II  load  current  (as  current  density)  (A  cm-2 ) 

Re  =  r£  +  rf  +  r“  electron  area  specific  resistance  of  the 
membrane  (£2  cm2) 

R0 2-  =  r^2  +  r®2_  +  r“2  oxygen  ion  area  specific  resistance 
of  the  membrane  (£2  cm2) 

Rh+  =  r^+  +  rf+  +  r“+  proton  area  specific  resistance  of  the 


X  =  RV 
Rl 

R'  = 


membrane  (£2  cm2) 
fKe  +R02-Re  +Rh+R02-  (£22 cm4) 

load  (given  as  area  specific  load)  (£2  cm2) 


(Re+Rl) 


net  electronic  area  specific  resistance  including 


internal  electronic  resistance  and  load  (£2  cm2) 

X'  =  Rh+  R'e  +  R0 2- R'e  +  Ru+  R0 2-  ( £22  cm4 ) 

Po2  oxygen  partial  pressure 

p h2  hydrogen  partial  pressure 

Ph2o  partial  pressure  of  H20 

R  ideal  gas  constant  (J  mol-1  K-1 ) 

kB  Boltzmann  constant  (J  K-1 ) 

T  absolute  temperature  (I<) 

F  Faraday  constant  (C  mol-1 ) 

z/  valence  of  species  i 

e  electronic  charge  (C) 

£02  Nernst  voltage  related  to  differences  in  /x02  across 
the  membrane  (V) 

Fh2  Nernst  voltage  related  to  differences  in  /zh2  across 
the  membrane  (V) 

Jo2  flux  of  02  through  the  membrane  (mol  s-1  cm-2 ) 

Jh2  flux  of  H2  through  the  membrane  (mol  s-1  cm-2 ) 

Jq2-  flux  of  oxygen  ions  through  the  membrane 

(mols-1  cm-2) 

JH+  flux  of  protons  through  the  membrane 

(mols-1  cm-2) 

Ljj  Onsager  coefficient  (mol2  cm-2  s-1  J-1 


Greek  letters 

oce  cathode/electrolyte  interface  electronic  conductiv¬ 

ity  (S  cm-1) 

Gq  anode/electrolyte  interface  electronic  conductivity 

(Son-1) 

of  electrolyte  electronic  conductivity  (S  cm-1 ) 

<7^2  cathode/electrolyte  interface  oxygen  ion  conductiv¬ 

ity  (S  cm-1) 

o“ 2  anode/electrolyte  interface  oxygen  ion  conductivity 
(Son-1) 

o'®2_  electrolyte  oxygen  ion  conductivity  (S  cm-1 ) 

<t^+  cathode/electrolyte  interface  proton  conductivity 

(Son-1) 

o^+  anode/electrolyte  interface  proton  conductivity 

(Son-1) 

0"®+  electrolyte  proton  conductivity  (S  cm-1 ) 

8C  cathode/electrolyte  interface  thickness  (cm) 

8a  anode/electrolyte  interface  thickness  (cm) 

fio2  chemical  potential  of  oxygen  (J  mol-1 ) 

/xh2  chemical  potential  of  hydrogen  (J  mol-1 ) 

/xh2o  chemical  potential  of  H20  (J  mol-1 ) 
fie  chemical  potential  of  electrons  (J  mol-1 ) 

0  electrostatic  potential  (V) 

jli  =  fii  +  Zie@  electrochemical  potential  of  species  i 
(J  mol-1 ) 

<p  =  -^  =  -^+0  electric  potential  (V) 


and 

X  +  me'  Xm~ 

(5) 

with  local  equilibria  given  by 

P'M  —  MMn+  +  =  A  Mn+  +  nAe 

(6) 

and 


fix  +  mjle  =  jxXm- 


(7) 


where  /xe  and  jle  are  respectively  chemical  and  electrochemical 
potentials  of  electrons.  The  /xe  in  an  aqueous  solution  is  a  well 
defined  quantity  even  if  the  actual  concentration  of  electrons  may 
be  negligible  so  that  negligible  current  is  carried  by  electrons. 
By  virtue  of  the  presence  of  mobile  cations  and  anions,  explicit 
incorporation  of  electron/hole  chemical  potentials  is  generally  not 
required.  Many  high  temperature  materials,  however,  exhibit  sub¬ 
stantial  mobility  for  typically  only  one  species.  For  example,  high 
temperature  cubic  phase  of  Zr02  stabilized  with  Y203  exhibits  sig¬ 
nificant  oxygen  ion  conduction  (via  a  vacancy  mechanism)  above 
about  400 °C.  Yttrium  (Y3+)  and  zirconium  (Zr4+)  ions  are  essen¬ 
tially  ‘frozen’  in  place  and  exhibit  no  detectable  conductivity  [3]. 
Description  of  local  equilibrium  in  these  materials  with  respect  to 
transporting  oxygen  ions  requires  the  incorporation  of  electronic 
defects  -  electrons  or  holes.  For  example,  in  an  oxygen  ion  conduc¬ 
tor  with  negligible  mobility  for  cations,  local  equilibrium  reaction 
may  be  given  by 

l02+2e'<s>02-  (8) 


with  local  equilibrium  given  by 

\^02  +  2Ae  =  Aq2-  (9) 


Local  equilibrium  between  oxygen  ions,  neutral  oxygen  and  elec¬ 
trons  (holes)  is  generally  assumed  regardless  how  low  might  be  the 
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electronic  conductivity,  which  is  still  much  greater  than  ionic  con¬ 
duction  of  Y3+  and  Zr4+.  There  are  a  number  of  solid  state  oxygen 
ion,  fluoride  ion,  chloride  ion,  sodium  ion,  lithium  ion,  potassium 
ion,  and  silver  ion  conductors  in  which  only  one  mobile  ion  exists. 
In  all  such  materials,  local  equilibrium  must  be  written  in  terms 
of  electrons  (holes).  This  is  an  important  difference  between  the 
extensively  investigated  aqueous  electrolyte  solutions  and  solid 
electrolytes,  and  yet  often  ignored  in  studies  on  solid  electrolytes 
[4]. 

There  are  many  alkaline  earth  perovskite  oxides  (e.g.  Y-doped 
BaCe03)  which  are  predominantly  oxygen  ion  and/or  electron  hole 
conductors  in  a  dry  atmosphere  [5].  When  heated  in  water  vapor, 
they  dissolve  water  thus  becoming  proton  conductors  [5].  Such 
materials  thus  can  transport  two  ionic  species  (H+  and  O2- )  in  addi¬ 
tion  to  electrons/holes.  Potential  applications  of  these  materials  are 
in  intermediate  temperature  (400-700  °C)  fuel  cells,  electrolyzers, 
and  hydrogen  separators. 

Transport  through  membranes  when  multiple  species  exhibit 
significant  mobilities  has  been  extensively  addressed  in  biologi¬ 
cal  systems.  In  such  systems,  significant  swelling  of  the  membrane 
occurs.  Many  reported  works  address  this  aspect  by  accounting 
for  various  frames  of  reference.  In  the  oxide  materials  considered 
here,  cation  immobility  implies  rigid  lattice  and  thus  a  fixed  frame 
of  reference.  Analysis  of  transport  in  biological  systems  has  been 
examined  using  what  is  known  as  ‘network  thermodynamics’  in 
which  analogies  between  electrical  circuits  and  transport  through 
membranes  are  drawn  [6].  Using  this  approach,  various  parame¬ 
ters  such  as  resistances  and  capacitances  can  be  given  in  terms  of 
transport  parameters  of  the  membrane  and  thermodynamics  of  the 
materials/systems. 

In  this  manuscript,  an  equivalent  circuit  approach  is  used  to 
examine  transport  through  mixed  proton,  oxygen  ion,  and  elec¬ 
tron  (hole)  conducting  membranes  under  steady  state  fuel  cell 
conditions.  Local  equilibrium  is  assumed  in  an  otherwise  non¬ 
equilibrium  state  [7].  Thus,  the  present  analysis  is  part  of  the 
general  topic  of  non-equilibrium  steady  state  [8,9].  By  virtue  of  the 
assumption  of  steady  state,  no  changes  in  ‘capacity’  occur.  Thus, 
the  equivalent  circuit  can  be  fully  described  in  terms  of  various 
resistances  and  internal  EMFs,  the  latter  related  to  the  chemical 
potentials  of  the  various  neutral  species.  This  allows  for  a  complete 
description  of  steady  state  in  terms  of  thermodynamic  potentials 
and  dissipative  components. 


2.  Analysis 

We  consider  a  solid  dense  membrane  (no  porosity  in  the  mem¬ 
brane)  capable  of  transporting  protons  (H+),  oxygen  ions  (O2-),  and 
electrons  (e)  or  electron  holes  (h).  Transport  of  electrons  or  holes 
will  be  treated  in  terms  of  a  single  electronic  conductivity  which 
embodies  transport  of  either  or  both  electronic  species.  Two  porous 
electrodes,  cathode  and  anode,  are  applied  on  the  two  surfaces  of 
the  membrane.  A  mixture  of  gaseous  02  and  H20  is  supplied  at 
the  cathode  and  a  mixture  of  gaseous  H2  and  H20  is  supplied  at 
the  anode.  An  external  load  is  connected  across  the  cathode  and 
the  anode.  Fig.  1  shows  a  schematic.  Cathode/membrane  (c)  and 
anode/membrane  (a)  interfaces  are  each  characterized  by  three  sets 
of  area  specific  resistances  (in  £2  cm2)— one  set  for  the  transfer  of 
H+  across  the  interfaces  (r£+,  r“+),  one  set  for  the  transfer  of  O2- 
across  the  interfaces  (rc9  ,  ra9  ),  and  one  set  for  the  transfer  of 
electrons/holes  across  the  interfaces  (r£,  r“) 2  Area  specific  resis¬ 
tances  for  ion  transfer  across  interfaces  can  be  given  in  terms  of 


2  The  interfacial  resistances  defined  here  do  not  include  extended  electrochemical 
zones  of  porous  electrodes. 


I  c  a  II 

J“H2  J“h2  Vh2  j“h2 


I - VVW - 1 

Load 

Fig.  1.  A  schematic  of  a  fuel  cell.  The  locations  of  the  chemical  potentials  of  H2 
just  inside  the  electrolyte  (/x ^ ,  /xj^ )  and  just  across  into  the  electrodes  (/xj^ ,  /x“2 ) 
are  labeled  in  the  schematic.  The  same  applies  to  the  other  chemical  potentials, 
electrochemical  potentials  and  (p. 


fundamental  charge  transfer  mechanisms  and  may  be  described 
phenomenologically  by  models  such  as  the  Butler-Volmer  model. 
Area  specific  resistances  for  electron  transport  across  interfaces 
(r£,  r“)  are  measures  of  direct  electron  transfer  between  the  elec¬ 
trodes  and  the  electrolyte.  For  many  situations,  electron  transport 
through  the  electrolyte  (and  across  interfaces)  is  negligible.  How¬ 
ever,  it  is  not  mathematically  zero.  This  means  rf ,  r“  may  be  very 
large,  but  are  not  infinite  [10].  Transport  through  the  membrane 
is  characterized  by  area  specific  resistances  for  H+  transport,  O2- 
transport,  and  electron/hole  transport  through  the  membrane. 
These  can  be  given  in  terms  of  the  respective  bulk  conductivities. 

Local  equilibrium  gives  for  all  regions  of  the  membrane  and  the 
adjacent  gas  phases  the  following  relations: 


h2  +  io2  ^  h20 

which  gives 
1 

Mh2  +  2^°2  =  /xh2o 
1h2  <^H+  +  e' 

which  gives 
1 

2 /^H2  —  Ah+  +  Ae 


(10) 


(11) 

(12) 


(13) 


and  Eqs.  (8)  and  (9)  describe  local  equilibrium  for  oxygen.  The 
above  (local)  equilibria  exist  everywhere  in  the  system,  including 
the  dense  membrane.  An  important  point  to  note  is  that  chemical 
potentials  of  gaseous  species  are  defined  in  fully  dense  membranes 
(no/negligible  porosity  in  the  bulk  of  the  membrane)  even  when 
the  amount  of  species  existing  as  gas  within  the  solid  is  negligible 
(or  mathematically  zero)  [11,12]. 

At  the  gas  phase/electrolyte  (membrane)  interface  I,  on  the  gas 
phase  side  of  the  interface,  local  equilibrium  gives 


M  h2  +  2^02  =  ^h20 


(14) 
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Let  I  be  the  cathode.  Then,  it  is  understood  that  for  a  fuel  cell,  the 
cathode  predominantly  has  O2  and  H20,  but  very  little  H2;  that  is 
Po2 ,  Ph2o  >  Ph2  •  At  the  gas  phase/electrolyte  (membrane)  interface 
II,  on  the  gas  phase  side  of  the  interface,  local  equilibrium  gives 

^2  +  \iA2  =  iA2o  n5) 

Let  II  be  the  anode.  Then,  it  is  understood  that  for  a  fuel  cell, 
the  anode  predominantly  has  H2  and  H20,  but  very  little  02;  that 
ispH2,PH2o  » Po2. 

Local  equilibrium  inside  the  electrolyte  (membrane),  just  under 


the  cathode/electrolyte  interface,  is  given  by 

Mh2  +  \r-o2  =  Mh2o  (16) 

Local  equilibrium  inside  the  electrolyte  (membrane),  just  under 

the  anode/electrolyte  interface,  is  given  by 

^h2  +  \^o2=^h2o  07) 

The  general  transport  equation  for  a  charged  species  is 
h  =  (18) 

Zit 

where  ot  is  the  conductivity  due  to  species  i  and  J;  is  the  current 
density  due  to  the  transport  of  species  i.  This  gives 

/h+  =  -^VAh+=-^V^h2-cth+V<p  (19) 

V-  =  %- v A02-  =  o2  -  ao2_  Vcp  (20) 

h  =  y  VAe  =  -aew<p  (21) 

where 


atomic  layers  thick)  which  is  different  from  both  the  solid  and  the 
gas  phases.  Indeed,  similar  approach  has  also  been  used  to  describe 
interfaces  between  liquid  and  vapor  using  non-equilibrium  molec¬ 
ular  dynamics  simulations  [14].  The  interface  defined  here  does  not 
include  relatively  long  range  effects  such  as  the  space  charge,  which 
can  be  readily  included  separately. 

Similarly,  for  the  anode/electrolyte  interface  (II),  the  area  specific 
resistances  are 


(24a) 


(24b) 


and 


(24c) 


where  8a  is  the  anode/electrolyte  interface  thickness,  and  cr“+ ,  o“2_ 
and  of  are  respectively  proton,  oxygen  ion  and  electronic  conduc¬ 
tivities  of  the  anode/electrolyte  interface.  In  general,  even  though 
the  interface  conductivities  and  thicknesses  are  not  separately  mea¬ 
surable,  interface  area  specific  resistances  can  be  measured,  at  least 
in  principle,  using  techniques  such  as  impedance  spectroscopy. 

Finally,  for  the  bulk  electrolyte  membrane,  the  area  specific 
resistances  are 


H+ 


(25a) 


and 


(25b) 


<P=  ~ 


Ae 

F 


(25c) 


An  important  point  to  note  is  that  what  one  experimentally  mea¬ 
sures  using  a  metal  wire  probe  is  cp  (electric  potential)  and  not 
the  electrostatic  potential,  0.  Much  of  the  literature  on  aqueous 
electrochemistry  only  refers  to  the  electrostatic  potential,  0  [13]. 

For  the  interface  at  I,  the  various  area  specific  resistances 
(Q  cm2)  are  as  follows: 


II 

+ 

(23a) 

rc  _ 
o2-  ac2 

O2- 

(23b) 

and 

rc  —  — 

(23c) 

where  8C  is  the  cathode/electrolyte  interface  thickness  (which  could 
be  in  the  nanometer  or  fractions  of  a  nanometer  range  and  the  inter¬ 
facial  region  need  not  be  physically  distinct),  and  a^+,  ocq2  and  of 
are  respectively  proton,  oxygen  ion  and  electronic  conductivities 
of  the  cathode/electrolyte  interface.  Although  the  exact  character 
and  the  thickness  of  the  interface  are  not  readily  subject  to  inde¬ 
pendent  experimental  determination,  some  important  features  are 
as  follows.  The  interface  region  should  include  one  (or  few)  atomic 
layers  of  the  solid  and  one  (or  few)  layers  of  the  gas  phase  (possibly 
adsorbed).  The  surface  layer  of  the  solid  must  always  be  ‘different’ 
from  the  bulk  region  due  to  the  presence  of  unsatisfied  bonds.  Thus, 
through  this  region,  the  transport  properties  must  also  in  general  be 
different  from  the  bulk.  Similarly,  the  adsorbed  layer  will  also  have 
different  properties  than  the  gas  phase.  In  this  discussion,  the  thick¬ 
ness  of  this  interface  is  meant  to  include  this  transition  region  (a  few 


where  i  is  the  membrane  thickness,  and  ofi. ,  oe\  and  of  are 
respectively  proton,  oxygen  ion  and  electronic  conductivities  of  the 
membrane.  These  in  principle  can  be  separately  measured  and  are 
typically  functions  of  composition  (defect  chemistry),  atmosphere 
and  temperature.  For  example,  oxygen  ion  conductivity,  o^2  ,  may 
be  given  in  terms  of  the  oxygen  vacancy  concentration  and  the 
oxygen  vacancy  diffusivity  (mobility).  An  example  is 


Dv ” 


1  _  4e2 _ o_  [1/" 

2-“4e  kBr 


(26) 


where  e  is  the  electronic  charge,  Dv»*  is  the  oxygen  vacancy  diffusiv- 

0 

ity  and  [\/"  ]  is  the  oxygen  vacancy  concentration.  Defect  chemistry 
relations  can  be  readily  given  for  specific  cases  for  all  species.  These 
are  material  and  atmosphere  dependent.  The  following  discussion 
is  applicable  regardless  of  the  details  of  defect  chemistry.  For  this 
reason,  no  particular  defect  mechanism  is  assumed  in  what  follows. 

Integration  and  rearrangement  of  Eqs.  (19)  through  (21)  (for  a 
one  dimensional  case)  assuming  transport  parameters  are  constant, 
independent  of  position,  gives 

o2  =  Mo2  +4F  {ro2-J02-  -  re'T 

=  Mo2  +  4Fr^2_/02-  -  4F  {*>'  -  <pc}  (27) 

l4\2  =  Mh2  _  2F  |r^+/H+  -  r£/e} 

=  m'H2  -  2Fr^+/H+  +  2F  [cpl  -  cpc }  (28) 

and 


<PC  =  <Pl  -  rceIe 


(29) 
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where  all  thermodynamic  potentials  (/xs  and  cp)  with  superscript  c 
refer  to  values  in  the  membrane  just  inside  the  cathode/electrolyte 
interface. 

Similarly,  for  the  anode/electrolyte  interface  (just  inside  the 
membrane) 

M&,  =  Mo,  -  4F  {ro2- V  -  re‘e} 

=  Mo,  -  4Fr“2_/02_  +4 F{<pa-  <p" }  (30) 


Mh2  =  Mh2  +  2 F  {r“+/H+  -  I'eF  } 

=  /_x[i2  +  2 FraH+IH+  -  2 F  {<pa  -  /  }  (31 ) 

and 

=  y>"  +  r“/e  (32) 

where  all  thermodynamic  potentials  (/xs  and  <p)  with  superscript  a 
refer  to  values  in  the  membrane  just  inside  the  anode/electrolyte 
interface. 

Eqs.  (27)  through  (32)  describe  thermodynamic  potentials 
within  the  electrolyte,  just  under  the  electrode/electrolyte  inter¬ 
faces.  The  most  important  feature  of  these  equations  is  that  they 
describe  thermodynamic  potentials  in  terms  of  ionic  and  electronic 
currents.  Thus,  by  virtue  of  the  presence  of  terms  containing  prod¬ 
ucts  of  electronic  resistances  and  electronic  currents,  that  is  of  the 
type  r“/e  and  r£/e,  it  is  to  be  noted  that  one  cannot  a-priori  assume 
electronic  resistances  as  being  infinite  (that  is  r“  oo  and  rce  ±  oo), 
as  is  commonly  done.  This  is  because  if  the  electronic  resistance 
is  infinite,  then  Ie  will  go  to  zero  resulting  in  indeterminate  terms 
(oo  x  0)  in  Eqs.  (27)  through  (32),  making  thermodynamic  poten¬ 
tials  indeterminate— and  thus  violating  local  equilibrium  [7,12]. 
Based  on  considerable  work  done  on  non-equilibrium  thermo¬ 
dynamics,  it  is  now  well  established  that  local  equilibrium  is 
applicable  in  virtually  all  practical  systems  of  interest,  such  as  the 
one  considered  here.  Thus,  local  equilibrium  must  be  valid— which 
here  means  electronic  resistance  cannot  be  assumed  to  be  infinite. 

2.1  The  equivalent  circuit 

An  equivalent  circuit  describing  steady  state  transport  can  be 
constructed  for  each  of  the  segments.  This  is  described  in  what 
follows. 


2.1.1.  The  cathode  I  electrolyte  interface  segment 

Rearrangement  of  Eq.  (28)  for  the  change  in  chemical  potential 
of  H2  across  the  cathode/electrolyte  interface  gives 


Mh2 


-ml 


2  F 


+  rfj+/H+  =  <pl  -  (pc  =  rceIe 


(33) 


Note  that  (pl  -  cpc  =  rcele  is  the  ‘voltage’  drop  or  ‘electric  potential’ 
drop  across  the  cathode/electrolyte  interface.  Write 


(a)  Cathode/Electrolyte 
Interface 


Mo2 

<pj 


v’rVVXA 

◄ —  ^o2~  Ea 

ww 

V)2- 


(b) 


Mo2 

<PC 


Electrolyte 

rAAAA-T? 


Ka/vh 


/„  2- 


WVb 

V)  2- 


Fet 


Mh2 

<P° 


(c)  Anode/Electrolyte 
Interface 


Mo: 


lo2~ 


A//- 


V)2- 


1«* 


LAA/V  hr 


F° 

hH, 


Fig.  2.  Equivalent  circuits  for  the  cathode/electrolyte  interface  (a),  the  electrolyte 
(b),  and  electrolyte/anode  interface  (c). 


Write 


co2 


>  0 


(37) 


Note  this  is  the  Nernst  potential  which  relates  the  change  in  free 
energy  of  transferring  one  mole  of  oxygen  from  /  to  c  (cathode). 
Thus, 


E02  +  ro2-  702-  =  re!e  =  <P'~<PC  (38) 

Fig.  2(a)  shows  the  equivalent  circuit  for  the  cathode/electrode 
interface  portion  of  the  cell.  Note  that 


cc  cc  ^H20  Mh20 
£02“£h2  = - 2 F - 


(39) 


In  general  the  magnitudes  of  E ^  and  Ec0 2  will  be  different.  The 
only  condition  for  which  they  would  be  the  same  will  be  if  /x^q  = 
/x^q,  which  would  be  a  very  unlikely  case.  In  general, /x^q  /xJ^q, 
even  when  the  partial  pressures  of  H20  may  be  the  same  at  the 
cathode  and  anode,  i.e.,  even  when  /Xh2o  =  ^h2o  (which  would  be 
the  case  if  fuel  and  oxidant  are  humidified  to  the  same  degree).  This 
aspect  will  be  discussed  later. 


E 


h2 


>  0 


(34) 


2.1.2.  The  electrolyte  segment 

Similar  analysis  as  above  for  the  electrolyte  segment  gives 


Note  this  is  the  Nernst  potential  which  relates  the  change  in  free  _  rei ^ 

energy  of  transferring  one  mole  of  hydrogen  from  c  to  /  (cathode). 

Thus,  For  hydrogen, 


(40) 


fh2  +  =  <p'-<pc 

Similarly, 


Mo,  -  Mo, 


+  rcQ2J02-  =(p'-<pc  =  rceh 


eie 


(35) 


Eel  - 
and 


Vh2-Vh2 
2  F 


>  0 


(36)  Ef  +  rH+^H+  =  rl'F  =  <pc  -  (pa 


(41) 


4F 


(42) 
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For  oxygen, 

2.2.  Nernst  voltages 

£e,  _^02-MS2 

b°2  ~  4 F  >U 

(43) 

The  Nernst  voltages  are  given  by 

and 

Eo2=E82+^2+£S2  =  g,n| 

60 

|  ^02  ^02 

(52) 

£02  +  ro2-I02-  =  reel/e  =  VC  -  Va 

(44) 

VP°J 

1  _  4  F 

Note  that 

and 

cel  cel  ^H20  “  ^H20 

(45) 

£h  2=E&2+Jtf2+ES2  =  §ln| 

j  ^h2  “  ^H2 

(53) 

E°2  £h2  =  2F 

^Ph2  ) 

‘  “  2  F 

In  general,  the  magnitudes  of  EjJ  and 

,  are  different.  Fig.  2(b) 

The  Eq2  and  Eh2  are  in  general  different.  They  will  be  equal  only 

shows  the  equivalent  circuit  for  the  electrolyte  portion  of  the  cell. 


2A.3.  Anode-electrolyte  segment 

Similar  analysis  for  the  electrolyte/anode  interface  gives 

cpa  -  (pn  =  r“Ie 


(46) 


if  Mh2o  =  ^h2o  (^at  same  degree  of  humidification  at  both 
electrodes). 

2.3.  Cell  ionic  and  electronic  area  specific  resistances  including 
membrane! electrode  interfaces 


For  hydrogen 

|Ltu  -/to 

Ea  _  '  H2  ^H2  q 

b»2  -  2 F  >U 

(47) 

and 

£H2  +  rH+^H+  —  r%Ie  =  <pa  - 

(48) 

For  oxygen 

Ea  _^02-^S2  >q 

b°2  -  4 F  >U 

(49) 

and 

£02  +  ra02-l02-  =  r“fe  =  <pa-  (Pn 

(50) 

Note  that 

pa  pa  MH20  -  M h2o 

eo2  Eh2  -  2  F 

(51) 

In  general,  the  magnitudes  of  2  and  ,  are  different.  Fig.  2(c) 
shows  the  equivalent  circuit  for  the  anode/electrolyte  interface 
portion  of  the  cell. 

Equivalent  circuit  for  the  complete  cell  at  open  circuit  is  merely 
a  series  connection  of  the  three  segments.  The  general  case  is  one 
where  an  external  load  is  connected.  This  is  shown  in  Fig.  3.  The 
various  parameters  of  the  equivalent  circuit  are  given  in  what  fol¬ 
lows. 


The  various  cell  area  specific  resistances  are 


!o2-=ro2-+r£-+ro2- 

(54) 

’H+  =  rH+  +  rH+  +  rH+ 

(55) 

■e  =  rce  +  rf  +  r“ 

(56) 

2.4.  The  current  densities 

The  various  current  densities  are  /H+.  I02 Je,  and  4  (4  is  the 
load  current,  given  as  current  density  per  unit  electrode  area).  In 
the  above  the  signs  of  the  various  current  densities  (with  respect 
to  Figs.  1  and  3,  and  with  the  positive  direction  through  the  mem¬ 
brane  from  left  to  right— from  the  cathode  through  the  membrane 
to  the  anode)  are  as  follows:  Iq2-  <  0,  /H+  <  0,  Je  >  0,  and  4  <  0.  The 
relationship  between  the  current  densities  is 

4  =  /02-+/h++4  (57) 

In  terms  of  the  various  resistances  and  current  densities,  the  cell 
voltage  is  given  by 

V  =  cpx  -  (pu  =  IeRe  =  -4^4  =  Eo2  +  I02-R02-  =  £h2  +  ih+eh+  (58) 

Analysis  of  the  equivalent  circuit  shows  that  the  cell  voltage, 
V=(pI-(pn  is  given  by 

V  =  tp]  -  tpn  =  E-  |/j.  |  Rcell  (59) 
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Fig.  3.  An  equivalent  circuit  for  a  fuel  cell  with  a  mixed  proton,  oxygen  ion  and  electron/hole  conducting  membrane.  The  block  arrows  show  the  locations  where  the  various 
potentials  within  the  electrolyte,  just  across  electrode/electrolyte  interfaces,  are  defined.  The  circuit  elements  between  the  filled  circles  at  the  ends  are  not  physically  separable. 
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where  \IL\  is  the  magnitude  of  the  load  current,  the  cell  open  circuit 
voltage  (OCV)  is 


E  = 


( 

V  R02-  ^H+ 


^02_^H+^e 

(l?H+^e  +  Rq 2- Re  +  ^H+^o2-  ) 


(60) 


and  the  cell  area  specific  resistance  is 


Rcew  — 


Rf^-R^+Re 

(^H+^e  +  R02-Re  +  ^H+^o2-  ) 


(61) 


For  convenience  of  a  shortened  notation  we  will  write  in  what 
follows 


X  =  RH+Re  +  RQ2-Re  +  ^h+^o2- 

The  short  circuit  current  density  is  given  by 


(62) 


k(s) 


(r^  +  tr) 

\R02-  rh+  j 


(63) 


Both  the  V  and  \IL  |  are  functions  of  the  load,  RL.  The  power  density 
at  any  load  is  simply  the  product  of  V  and  \IL\. 

The  OCV  (Eq.  (60))  is  effectively  the  cell  potential  as  defined 
in  the  Goldman-Hodgkin-Katz  equation  [1,2],  which  is  the  cell 
potential  (an  experimentally  measured  quantity)  in  terms  of  the 
individual  Nernst  voltages  and  transport  parameters.  In  much  of  the 
literature  on  solid  electrolytes,  the  cases  on  interest  involve  a  single 
Nernst  voltage.  An  example  of  a  system  with  two  Nernst  potentials 
is  described  here  for  a  solid  electrolyte  membrane  transporting  two 
ions.  In  biological/aqueous  systems,  the  GHK  equation  is  given  in 
terms  of  the  properties  of  the  membrane  and  ignoring  the  effects 
of  interfaces.  If  the  membrane  is  thin,  however,  interface  contribu¬ 
tions  dominate  transport  resistance  through  the  membrane.  Thus, 
for  thin  membranes,  it  is  necessary  to  include  interface  effects.  Eqs. 
(60)  and  (61)  given  here  do  take  into  account  interface  effects  and 
thus  are  applicable  to  membranes  of  any  thickness. 

For  a  predominantly  O2-  conductor  with  negligible  FI+  con¬ 
duction  and  negligible  electron  transport,  Eq.  (59)  reduces  to  V  = 
(pl  -  (pn  =  Eo2  —  |/L  |  Rq2~  .  For  a  predominantly  FT  conductor  with 
negligible  O2-  conduction  and  negligible  electron  conduction,  Eq. 
(59)  reduces  to  V  =  (pl  -  cpn  =  Fh2  -\h  \  Rh+-  These  equations  are 
commonly  used  when  describing  a  solid  oxide  fuel  cell  (SOFC)  based 
on  an  oxygen  ion  conductor  (e.g.  YSZ)  [15]  or  a  proton  exchange 
membrane  fuel  cell  (PEMFC)  using  a  proton  conductor.  Eqs.  (58) 
through  (63)  and  the  equivalent  circuit  describe  the  behavior  of 
mixed  proton,  oxygen  ion,  electron  (hole)  conducting  membranes, 
such  as  various  doped  alkaline  earth  perovskites  (e.g.  Y-doped 
BaCe03  or  Y-doped  BaZr03).  These  equations  also  take  into  account 
transport  of  all  three  species  (H+,  O2-  and  electrons/holes)  across 
electrode/membrane  interfaces. 


2.5.  Transport  under  open  circuit  condition  (OCV) 

At  OCV,  the  load  current  is  zero.  When  the  load  current  is  zero 
(141  =  0),  that  is  when  the  load  resistance  is  infinite,  no  net  current 
flows  through  the  cell.  In  such  a  case,  Iq2-  +  IH+  +  Ie  =  0.  Solution  to 
the  transport  equations  for  proton  and  oxygen  ion  current  densities 


gives, 

Re  +  Rh+  Re 

lo2~  =  X  L°2  +  X tH2 

(64) 

and 

Re  Re+  Rq  2- 

jH+  =  X  h°2  X  tH2 

(65) 

Eqs.  (64)  and  (65)  satisfy  Onsager  reciprocity  relations,  as  is  to 
be  expected.  The  corresponding  electronic  current  density  is  given 


by3 

j  _  RH+  p  ,  V-  r 

Ie  -  ~1CEo2  +  —  Eh2 


(66) 


The  ion  fluxes  (in  mol  cm-2  s-1 )  are  related  to  the  current  den¬ 
sities  as  follows: 


t02-  =  -2  FJq2-  (67) 

and 


Jh+=FJh+  (68) 

There  is  no  net  current  flowing  through  the  cell.  The  ion  fluxes 
are  equivalent  to  the  corresponding  neutral  species  ‘permeation’ 
fluxes.  Thus, 

Iq2-  =  — 2F/02-  =  -4FJq2  (69) 

and 


Ih+=FJh+=2FJH2  (70) 

where  Jo2  and  Ju2  are  permeation  fluxes  of  O2  and  H2  respectively 
across  the  membrane  (cell). 

Eqs.  (69)  and  (70)  thus  can  be  rearranged  to  read 


Re  +  Fh+  Re  a 

°2  -  16 F2X  ^°2  8 F2XA/XH2 

(71) 

and 

Re  a  Re  T  V-  a 

Jli2  ~  8F2X  AM°2  4F2X  A/lH2 

(72) 

where 

A^o2  =  M o2  -  Mo2 

(73) 

and 

A^h2  =  M h2  -  M h2 

(74) 

Note  Eqs.  (71)  and  (72)  are  of  the  form 

Jo2  =  -Lo2o2Afio2  -Lo2h2A^h2 

(75) 

J h2  =  -Th2o2A/xo2  -Lh2h2^/xh2 

(76) 

where 

Re  +  FH+ 

°2°2  “  16F2X 

(77) 

Re  +  R02- 
Lh2H2=  4F2X 

(78) 

Fh202  =  Fo2h2  =  8F2X 

(79) 

Note  that  Eqs.  (71),  (72),  (75),  and  (76)  are  the  Onsager  equations 
and  the  cross  terms  satisfy  the  Onsager  reciprocity  relations.  In 
these  equations,  the  driving  forces  are  given  not  as  gradients  but 
as  differences  in  chemical  potentials  across  the  membrane,  that  is, 
A fio2  =  /Xq2  -  /Xq2  and  A/zh2  =  /x}^  -  /xj^.  Thus,  the  Ly-s  defined 
here  contain  an  additional  length  dimension  in  the  denominator 
compared  to  the  usual  definition  [6].  Eq.  (79)  shows  that  the  reci¬ 
procity  relations  are  obeyed  not  only  for  the  bulk  membrane  but 
also  include  interface  transport  properties  embodied  in  Rq2-  (i.e. 
ro2->  ro2-)’^H+  (i*e.r£+,  r“+)andFe(i.e.r£,  r“).  It  was  assumed  in  the 
foregoing  that  transport  properties  of  the  membrane  are  position- 
independent.  If  the  transport  properties  do  depend  on  position,  the 


3  It  is  easy  to  show  that  if  the  electrolyte  material  is  a  predominantly  O2-  con¬ 
ductor  with  negligible  H+  conduction,  then  Iq2-  and  Ie  are  coupled.  On  the  other 
hand  if  the  electrolyte  material  is  a  predominantly  H+  conductor  with  negligible 
O2-  conduction,  then  JH+  and  Ie  are  coupled. 
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cell  membrane  may  be  discretized  into  as  many  segments  as  needed 
for  any  desired  level  of  accuracy  [12].  As  long  as  steady  state  is 
maintained,  equations  of  the  form  given  in  (75)  through  (79)  con¬ 
tinue  to  be  valid.  Also,  interface  parameters  are  naturally  included 
in  the  transport  equations.  Thus,  the  generalized  form  of  Onsager 
transport  equations  as  given  in  terms  of  A/xs  across  the  membrane 
remains  valid,  inclusive  of  interface  effects  even  when  transport 
parameters  depend  upon  local  composition  (and  thus  on  position). 


2.6.  Chemical  potentials  in  the  electrolyte  just  under 
electrode/electrolyte  interfaces 


Substitution  of  the  equations  for  current  densities  into  Eqs.  (27), 
(28),  (30),  and  (31 )  gives  the  chemical  potentials  /x^  ,  /xf^ ,  /x^ ,  and 
/x^2  as  follows: 


^02  =  Mo2  +  4F 

f  [r^2_(Re  +  Rh- )  +  r‘eRH - ]  _  [r‘ Re  - rp?Q2_]  'j 

x  { - x - E°2+ - x - E"2j 

(80) 

Mh2  =  Mh2  -  2F 

f[rcH+Re-r^RH+]  [i*  (Re  +  Rq2-  )  +  rceR02-  ]  } 

x  { - x - E°2 - x - E"2) 

(81) 


Mo2  =  Mo2  -  4 F 

r  [r“2  (Re  +  R^)  +  n!R„-\  [r“  Re-r?R0 2_]  'j 

x  \ - X - E°2+ - X - Eh2j 

(82) 

I^h2  =/4+2F 

f  [ra  Re  -  r“RH-  ]  [ra  (Re  +  R0 2- )  +  raeRQ2  }  } 

x  | - x - E°2 - X - Eh2/ 

(83) 

Eqs.  (80)  through  (83)  describe  chemical  potentials  of  02  and  H2 
in  the  electrolyte,  just  inside  the  electrode/electrolyte  interfaces  at 
OCV  in  terms  of  parameters  which  are  measureable  experimentally, 
at  least  in  principle.  An  important  point  to  note  is  that  in  general 
the  terms  in  *{}’  are  not  zero.  Thus,  the  chemical  potentials  must 
exhibit  abrupt  changes  across  interfaces  [12,16].  That  is,  in  general, 

^H2  ^  ^H2’  ^H2  ^  ^!v  ^02  ^  ^02  anC*  ^02  ^  ^02* 

In  many  reported  works  on  solid  electrolytes,  it  has  been 
commonly  assumed  that  chemical  potentials  equilibrate  across 
interfaces.  The  present  work  shows  that  such  equilibration  is  not 
to  be  expected.  In  fact,  as  shown  previously  [16],  equilibration  of 
chemical  potential  across  interfaces  is  generally  not  valid  even  if 
the  electrode  is  reversible.  That  is,  reversible  electrode  does  not 
mean  equilibration  of  the  chemical  potential  of  neutral  species 
across  electrode/electrolyte  interfaces  as  is  often  assumed,  when 
there  is  a  finite,  nonzero  chemical  potential  difference  across  the 
membrane.  Indeed,  recent  work  on  other  systems  (vapor/liquid 
interfaces)  using  non-equilibrium  molecular  dynamics  simulations 
has  also  shown  that  there  are  abrupt  chemical  potential  jumps 
across  interfaces  [  17].  It  is  easily  seen  that  the  cp  also  exhibits  abrupt 
changes  across  the  two  interfaces.  Fig.  4  shows  a  schematic  of  the 
variation  of  /xh2  ,  /xo2  and  (p  through  the  membrane  and  discontin¬ 
uous  changes  across  interfaces. 

Abrupt  changes  in  chemical  potentials  and  electric  potential 
across  interfaces  represent  a  manifestation  of  the  existence  of 
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Fig.  4.  Schematic  variations  of  /xh2  ,  H o2 .  and  (P-  Discontinuous  changes  in  chemical 
potentials  and  electric  potential,  <p,  occur  across  the  interfaces  in  all  functions. 


rather  thin  interfacial  regions.  This  means  that  chemical  potentials 
and  electric  potential  change  smoothly  within  the  individual  phases 
(gas  and  solid),  but  change  abruptly  at  the  interface.  Within  the 
interfacial  region  itself,  it  is  expected  that  chemical  potentials  and 
electric  potential  vary  continuously.  The  thinner  is  the  interfacial 
region,  the  more  sharply  do  the  chemical  potentials  and  the  electric 
potential  vary  within  the  interfacial  region.  Thus,  on  a  microscopic 
scale  (micron  or  a  fraction  of  micron),  discontinuous  changes  in 
chemical  potentials  and  electric  potential  occur  across  interfaces. 

Substitution  for  /Xq2  and  /xj^  in  Eq.  (16),  and  /Xq2  and  /xf^  in 
Eq.  (17)  gives  the  corresponding  chemical  potentials  of  H20  in  the 
electrolyte,  just  under  the  electrode/electrolyte  interfaces,  namely, 


^h2o  ~  ^h2o  - 


lrH+Re  +  rcQ2_(Re  +  RH+)] 


Eo2 


[ru+(Re  +  Rq2-  )  +  Fq2- Ee 


-£h. 


(84) 


^h2o  -  ^h2o  + 


l^+Re  +  raQ2_(Re  +  RH+)] 


E0? 


[r^+  ( Re  +  Rq  2- )  +  r“  2 


-Eh2 


(85) 


Eqs.  (84)  and  (85)  show  that  the  chemical  potential  of  H20  must  also 
exhibit  abrupt  changes  across  interfaces  since  in  general  terms  in 
*{}’  are  nonzero.  This  result  also  shows  further  interesting  features. 
Suppose  both  the  cathode  and  the  anode  are  humidified  to  the  same 
degree,  that  is  suppose  =  /x^q.  Even  then,  hlMiQ  +  /x“  i0.  This 
means  there  will  be  a  gradient  in  /x h2o  through  the  membrane  even 
when  the  same  partial  pressure  of  H20  exists  in  the  two  gas  phases 
(cathode  and  anode). 


2.7.  Currents  under  load 

When  an  external  load  is  connected  so  that  a  finite  current  (4) 
flows  through  the  external  circuit,  oxygen  ion  and  proton  current 
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densities  are  given  by 

Re  +  Rh+  Re  F  Rii+Re  f 

v-  = - x — E° 2  +  XEh2  +  L 

and 


/H+  =  JL En ,  - 


Re  +  Rn2~ 


02 


X 


-Eh2  + 


R02-Re 

X 


(86) 


(87) 


and 


L°2h2  Lh2°2  8F2(Rh+  +  Rq2- ) 


(101) 


In  this  case,  the  origin  of  coupling  is  purely  ionic.  In  such  a  case 
^o2o2^h2h2  _ (^h2o2)  =  °  (102) 


Eqs.  (86)  and  (87)  can  also  be  written  in  terms  of  the  load  resistance, 
Rl,  as 


R'e  +  Rh+  R'e 

= - x' — E°2  +  X7^2 


"o- 

and 


r  ReF  Re+Rn 

lH+  =  ^h02  ~ 


x' Eh2 

where  we  have  written 

ReRL 


K  = 


Re  +  RL 


(88) 


(89) 


(90) 


which  is  the  net  electronic  area  specific  resistance,  including  the 
external  load.  Also, 


Thus,  the  degree  of  coupling  between  H2  and  02  transport  through 
the  membrane  is  a  function  of  the  net  electronic  resistance  (includ¬ 
ing  internal  and  external)  and  varies  between  no  coupling  for  zero 
electronic  resistance  to  maximum  coupling  for  infinite  electronic 
resistance. 

The  chemical  potentials  of  H2,  02,  and  H20  in  the  electrolyte, 
just  under  the  electrode/electrolyte  interfaces,  are  the  same  as 
(80)-(85)  with  Re  replaced  by  R'e  which  includes  the  load  resistance, 
Rl  (thus  also  replacing  X  by  X'). 

2.8.  Significance  of  equations  describing  chemical  potentials 
inside  the  electrolyte  as  a  function  of  various  current  densities  and 
transport  parameters 


X  — Ry]+ R'e Rq2~  R'e  R}\+ Rq2-  (91) 

Note  that  Eqs.  (88)  and  (89)  are  similar  to  Eqs.  (64)  and  (65)  and 
satisfy  the  Onsager  reciprocity  relations.  Fluxes  of  02  and  H2  are  also 
similar  to  Eqs.  (71 )  and  (72)  with  Re  replaced  by  R'e,  and  X  replaced 


by  X',  and  are  given  by 

r  R'e  T  ^h+  a  R'e  a 

°2  ~  16F2X'  /X°2  8 F2X' A/iH2 

(92) 

and 

1  R'e  a  ^ e  T  Ro2~  a 

J n2  "  8F2X'  A/X°2  4F2X'  A/Xh2 

(93) 

Eqs.  (92)  and  (93)  can  still  be  considered  as  ‘permeation’  fluxes 
wherein  part  of  the  coupling  electron  flux  is  through  the  external 
load.  That  is,  Eqs.  (92)  and  (93)  represent  permeation  fluxes  for 
the  combined  system  comprising  the  cell  and  the  load.  The  corre¬ 
sponding  transport  coefficients  are  the  same  as  Eqs.  (77)  through 
(79)  with  Re  replaced  by  R'e  and  X  replaced  by  X',  namely 


T  f  Re+  ^H+ 

°2°2  16F2X' 

(94) 

_  R'e  +  V- 

H2H2  _  4F2X' 

(95) 

Lh2o2  —  lo2h2  —  SF2X' 

(96) 

Finally,  it  is  easily  seen  that 

lo2o2lh2h2  -  (lh2o2)2  ^  0 

(97) 

as  required  [18]. 

Various  limiting  cases  of  the  above  equations  can  be  examined, 
all  satisfying  Onsager  reciprocity  equations.  In  the  short  circuit 
limit,  R'e  =  0,  for  which 


i '  —i'  —  n 

LH202  -  l02H2  -  u 


(98) 


That  is,  there  is  no  coupling  between  H2  and  02  fluxes  at  short 
circuit.  Consider  now  R'e  -*  oo,  i.e.  no  load  applied  and  negligible 
electronic  conduction  through  the  membrane.  Then  the  corre¬ 
sponding  transport  parameters  are 


Lo2o2  16F2(Rh+  +  Rq2-  ) 
Lh2H2  =  4F2(Rh+  +  Rq2_) 


(99) 

(100) 


Eqs.  (80)  through  (83)  describe  chemical  potentials  of  02  and 
H2  in  terms  of  the  various  current  densities  (/02-,  /H+»  and  4)  and 
transport  parameters  at  OCV.  The  same  equations  are  applicable 
under  load  with  Re  replaced  by  R'e  and  X  replaced  by  X'.  A  key 
point  to  note  is  that  recognition  and  incorporation  of  local  equi¬ 
librium  has  made  it  possible  to  determine  chemical  potentials  of 
neutral  species  within  the  electrolyte.  These  equations  show  that 
for  a  single  fuel  cell,  since  the  ionic  current  densities  (/Q 2-  and 
/H+ )  are  in  the  opposite  direction  to  the  electronic  current  density 
(Je),  the  /zh2  and  /z02  in  the  electrolyte  are  always  mathematically 
bounded  by  values  at  the  electrodes  (that  is  by  ^  and  /z}^ ;  and 
by  /Zq2  and  /jl1*  )  [12].  However,  when  many  cells  are  connected 
in  series,  in  an  unbalanced  cell  the  ionic  and  electronic  currents 
can  be  in  the  same  direction  (the  electronic  current  changes  direc¬ 
tion).  In  such  a  case,  the  chemical  potentials  of  H2  and  02  within 
the  electrolyte  are  not  mathematically  bounded.  There  are  some 
situations  in  which  stack  failure  can  occur  as  described  previously 
[16,19]. 

2.9.  Possible  experimental  procedure  for  the  verification  of  the 
transport  equations 

The  analysis  given  here  assumes  that  the  transport  parameters 
(Rq2-  ,  Rh+  and  Re)  are  independent  of  gas  phase  compositions  and 
of  the  net  current  transporting  through  the  membrane.  In  many 
experimental  situations,  this  may  not  be  the  case.  Nevertheless,  if 
this  assumption  is  approximately  valid,  an  experimental  approach 
can  be  envisioned  for  the  measurement  of  transport  parameters.  In 
what  follows,  this  is  illustrated  using  Eqs.  (92)  and  (93). 

The  approach  consists  of  applying  a  known  load  (ftL)  across  the 
cell.  The  cathode  is  exposed  to  a  mixture  of  02  and  H20.  The  anode 
is  exposed  to  a  mixture  of  H2  and  H20.  The  corresponding  chem¬ 
ical  potential  differences  of  02  and  H2,  namely  A/z02  and  A/zh2, 
can  be  determined  based  on  the  compositions  (partial  pressures) 
of  the  gases.  By  conducting  chemical  analysis  of  incoming  and  exit¬ 
ing  gases  at  the  two  electrodes  and  measuring  gas  flow  rates,  the 
net  fluxes  of  02  and  H2,  namely  Jo2  and  Jh2>  can  be  experimentally 
measured  as  a  function  of  the  imposed  A/zo2  and  A/zh2.  Then,  the 
gas  mixtures  of  the  two  electrodes  are  changed  such  that  only  one  of 
the  chemical  potential  differences  is  altered  at  a  time.  For  example, 
consider  changing  A/z02  but  keeping  A/zh2  fixed.  The  net  fluxes  of 
02  and  H2,  namely  Jo2  and  Ju2,  are  again  measured.  This  procedure 
is  repeated  for  various  values  of  A/zo2  for  a  fixed  value  of  A/zh2. 
The  rate  of  change  Jo2  as  a  function  of  A/zq2  for  a  fixed  A/zh2  gives 
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Lq2q2  (Eq.  (94)).  Similarly,  the  rate  of  change  of  JHl  as  a  function  of 
A iiq2  for  a  fixed  A/zh2  gives 

A  similar  series  of  experiments  is  conducted  wherein  the  A/xh2 
is  systematically  varied  while  fixing  A/x02.  In  this  case,  the  rate  of 
change  of  Jo2  as  a  function  of  A/xh2  for  a  fixed  A/xo2  gives 
and  the  rate  of  change  of  Ju2  as  a  function  of  A/xh2  for  a  fixed  A/x02 
gives  H  (Eq.  (95)).  The  L'H  0  and  L'0  H  determined  from  the 
two  sets  of  experiments  should  be  equal  to  each  other  and  given  by 
Eq.  (96).  The  above  describes  a  possible  experimental  approach  for 
verification  and  general  application  of  the  equations  derived  here, 
and  for  the  verification  of  Onsager  reciprocity  relations  for  transport 
through  mixed  O2-,  H+,  and  electron  (hole)  conductors.  From  these 
measurements,  transport  parameters  for  the  membrane,  namely, 
Rq2~,Rh+,  and  Re  can  be  determined. 

3.  Summary 

The  present  manuscript  has  examined  steady  state  transport 
through  mixed  proton,  oxygen  ion,  and  electron  (hole)  conductors. 
Transport  through  membranes  in  a  fuel  cell  mode  was  exam¬ 
ined.  Transport  properties  of  all  transporting  species  (H+,  O2-,  e 
or  h)  across  cathode/electrolyte  and  anode/electrolyte  interfaces 
were  incorporated  into  the  analysis  and  a  simple  equivalent  circuit 
was  presented.  Only  steady  state  was  examined.  This  allowed  for 
the  use  of  an  equivalent  circuit  in  terms  of  internal  EMFs  (mea¬ 
sures  of  chemical  potential  differences)  and  resistors  (dissipative 
components),  without  the  need  for  elements  describing  ‘capac¬ 
ity’.  Ion  transport  properties  across  interfaces  are  measures  of 
electrochemical  charge  transfer  reactions  at  the  physically  sharp 
electrode/electrolyte  interfaces  (not  including  the  extended  porous 
electrode).  Electron  transport  properties  across  interfaces  are  mea¬ 
sures  of  direct  electron  transfer  between  the  electrodes  and  the 
electrolyte.  The  open  circuit  voltage  of  cell/membrane  is  given  in 
terms  of  chemical  potentials  of  H2  and  02  at  the  two  electrodes 
and  transport  properties  of  the  membrane  including  interfaces. 
The  resulting  equation  is  similar  to  the  Goldman-Hodgkin-Katz 
equation  used  in  cell  physiology  describing  ion  transport  across 
cell  walls.  Transport  through  the  membrane,  inclusive  of  interface 
transport  properties,  obeys  Onsager  relations.  In  many  situations, 
membranes  can  be  very  thin.  In  such  cases,  interface  effects  domi¬ 
nate  transport.  The  present  analysis  shows  that  Onsager  reciprocity 
relations  are  valid  with  interface  parameters  included  in  the  trans¬ 
port  equations.  In  the  analysis,  transport  properties  through  the 


bulk  membrane  were  assumed  to  be  position-independent.  If  they 
are  position-dependent,  the  entire  cell  can  be  discretized  into 
as  many  elements  as  required  for  a  specified  accuracy  [12].  The 
Onsager  equations,  also  accounting  for  interface  transfer,  continue 
to  be  valid.  The  present  manuscript  has  examined  transport  of  two 
independent  species  (H  and  O).  This  resulted  in  two  Nernst  poten¬ 
tials  and  three  parallel  circuit  elements  (including  one  for  electronic 
path)  for  any  given  segment  (three  segments— two  interfacial  and 
one  bulk)  in  the  equivalent  circuit.  The  generalization  of  this  result 
is  that  if  there  are  N  independent  transporting  ionic  species,  there 
will  be  N  Nernst  potentials  and  (N+ 1 )  parallel  circuit  elements  for 
any  given  segment. 
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